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I. INTRODUCTION 


Abstract— In this research, the breaker index equations are formulated 
using the kinematic free surface boundary condition. By substituting the 
potential velocity equation for the solution of Laplace's equation in this 
equation, it is obtained the wave amplitude function equation. From the 
wave amplitude function equation two breaker indices are extracted, 
they are the breaker length index which is the ratio between the breaker 
height and the breaker length; and the breaker depth index which is the 
ratio between the breaker height and the breaker depth. The next 
breaker index, which is a ratio between breaker depth and breaker 
length, is obtained from the wave number conservation law. Consistency 
testing of the three breaker index equations obtained shows that there is 
consistency in the three equations. Consistency testing is done by using 
the connectivity equation, where a breaker index is the product of the 
multiplication of the other two breaker indexes. The breaker height 
index, which is the ratio between the breaker height and the deep water 
wave height, is obtained by substituting the breaker length index in the 
energy conservation equation. Thus the breaker height equation is 
obtained which works a function of the breaker length at the breaking 
point. With the availability of the four breaker indexes, the breaking 
parameter can be calculated easily. 


breaker index is hereinafter referred to as the breaker 


The breaking parameter is the wave characteristics at the 
time of breaking. There are three breaking parameters, 
namely the breaking wave height or also known as the 
breaker height, the water depth at the breaking point 
which is called the breaker depth and the wavelength at the 
breaking point which is called the breaker length. Breaker 
index is a ratio between breaking parameters. The breaker 
depth index is the ratio between the breaker height and the 
breaker depth, the breaker length index is the ratio between 
the breaker height and the breaker length. There is a 
breaker index which is rarely researched or used, namely 
the ratio between breaker depth and breaker length; this 
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depth-length index. Breaker height index is the ratio 
between breaking wave height and deep water wave 
height. 

Several explicit equations to the breaker depth index were 
put forward by, among others, McCowan (1894), Weggel 
(1972), Galvin (1969), Collins and Weir (1969), Madsen 
(1976) and Smith and Krauss (1989). 

There is a critical wave steepness criterion, which is the 
wave steepness before breaking which is the ratio between 
the wave height and wavelength before breaking. This 
criterion can also be considered as a breaker length index. 
Two researchers who put forward the criterion of critical 
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wave steepness are Michell (1894) and Toffoli et al 
(2010). The breaker length index is in the form of an 
implicit equation which is a function of another breaker 
index, namely the Mieche equation (1944). This equation 
connects the breaker length index with the breaker depth- 
length index. Many researchers have developed the 
breaker length index by modifying the Mieche equation 
(1944), including Battjes and Jansen (1978), Ostendorf and 
Madsen (1979) and Rattanapittikon and Shibayama (2000). 
The various breaker index equations are formulated 
separately, where the breaker depth index is formulated 
only by examining the breaker depth index, the breaker 
length index is formulated only by examining the breaker 
length index, as well as the breaker height. There should 
be a relationship or connectivity between breaker indexes, 
where the value of one breaker index is related to the value 
of another breaker index. In the breaker index equations 
obtained in this study, an examination of the linkages 
between breaker indexes was carried out using the 
connectivity equation. As a result, it is found that there is 
connectivity between breaker indexes. 


II. THE VELOCITY POTENTIAL EQUATION 
The complete velocity potential solution of the Laplace 
equation (Dean (1991)) using the variable separation 
method is: 
o(x,z,t) = G (cos kx + sin kx) coshk(h + z) sin ot 

(1) 
Where x is the horizontal axis, z is vertical axis and t is 
time, while h is water depth. G, k dan ø is called the wave 
constant, G is the energy transmission rate, k is wave 


21 : 2m . 
number, k = Fa where L is a wavelength and o = 7 ÍS 


angular frequency while T is for wave period. 


The three wave constants need to be determined by their 
value or equation. As shown in (1), the velocity potential 
consists of two components, namely the cos and the sin 
component. In both components there is a point 
where cos kx = sinkx, these points are called 
characteristic points. Analysis of wave constants will be 
easier if done at characteristic points, where the constants 
obtained will satisfy both wave components. By using only 
components cos kx then the velocity potential equation 
becomes, 

o(x,z,t) = G cos kx cosh k(h + z) sin ot ...(2) 


Note that the G in (2) has a double value. 
HI. THE WAVE AMPLITUDE FUNCTION 
EQUATION 


The wave amplitude function is the relationship between 
wave amplitude and other wave constants. This equation is 
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obtained by substituting the velocity potential for the 
kinematic free surface boundary condition and the 
equation obtained is integrated with time f to obtain the 
elevation water surface equation. 


As a kinematic free surface boundary condition, a 
weighted kinematic free surface boundary condition from 


Hutahaean (2022) is used in the form, 


an _ = dn 
Y2 ie = Wn Uns, ...(3) 


Y2 is a coefficient that is greater than 1, but the value of 
this coefficient has no effect on the breaker index. n(x, t) 
is the water surface elevation relative to the still water 
level, w, surface vertical water particle velocity u, surface 
horizontal water particle velocity. 

By using the velocity potential equation, 


= Gk sin kx coshk(h + z) sin øt 
Uy = Gk sin kx cosh k(h + n) sin ot 


d 
w(x,z,t) = a 


= —Gk cos kx sinh k(h + z) sin ot 
W, = —Gk cos kx sinh k(h + 7) sin ot 


Substitute these velocity equations into (3) and work on 


the characteristic points, 


dn Gk 
— = —— coshk(h +n) 
Y2 


(tanh k(h+n) + an) coskxsinot ....... (4) 


As a periodic function, then 
d 
Gk coshk(h +n) (tanh k(h+n) + =z) 


= constant 
Then integration (4) can be completed by integrating the 
elements of sin ot, and is obtained water surface elevation 
equation of, 


Gk 
n(x,t) = — coshk(h+n) 
OY2 


(tanh k(h+n) + =) cos kx cosot 


Wave amplitude is defined as, 
Gk dn 
A = — coshk(h+n) (tanh k(h+n)+ =) 
OY2 dx 


Considering that G has a double value, then 


dn 
ay; coshk(h +n) (tanh k(h+n) + N 


A is the wave amplitude, the water surface elevation 
equation becomes, 
n(x, t) = A cos kx sin ot 
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At a characteristic point of space and time, 7 = f dan = = 
— a wave amplitude function becomes 
z nk (n+) (t nk (n+5) | 
= Jay, cos 5 an 3 7 
re (5) 


In the wave amplitude function there is an element of 


tanh k (h + $) which has a constant value. The constant 


value of the function of tanh k (h + $) is reached at a 


large value of h, ath = họ, ho is the deep water depth 
where the influence of waves still reaches the bottom of 
the waters even though it is very small or close to zero. 
The constant value is, 
Ao 

tanh ko (ho + =) =1 
Therefore, in deep waters even at depths greater than ho, 
the waves seem to move on water depth hy, where 

Ao 
ko (ho + 3) = On 
0 is the deep water coefficient whose value needs to be 
determined. However, based on the wave number 
conservation law discussed in section (6), this equation 
applies at all depths, so it can be written in general as: 


A 
k (h + $) -v (6) 
The determination of the value @ is to use tanh Or = 1. 


However, this criterion is very dependent on the level of 
accuracy used, as shown in Table (1). 


Table (1) Value tanh 67 


tanh 67 

0 (5) (6) (7) 
1.95 0.99999 0.999990 0.9999905 

2 0.99999 0.999993 0.9999930 
2.05 0.99999 0.999995 0.9999949 
2.1 1 0.999996 0.9999963 
2.4 1 0.999999 0.9999994 
2.45 1 1.000000 0.9999996 
2.75 1 1.000000 0.9999999 
2.8 1 1.000000 1.0000000 


In Table (1), (5), (6) and (7) is the number of digits after 

the decimal point. 

- With an accuracy of 5 decimal places, it is found that 
tanh Oz = 1 at 0 = 2.1 

- With an accuracy of 6 decimal places, it is found that 
tanh 0r = 1 at 0 = 2.45 

- | With an accuracy of 7 decimal places, it is found that 
tanh 0r = 1 at 0 = 2.80 
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Due to this condition, the value of 0 cannot be determined 
only by using the criterion of tanh Or = 1, thus additional 
criteria are needed. The more precise determination of 
value 0 will be done in section 9. 


IV. Equation of Breaker Length Index r, 
b 


In (5) there is a breaking characteristic, when 
A kA 
tanhk (h+) -#50 00 2 (7) 


In this condition the wave amplitude function is zero. In 
(7) substituted (6) in the 1* term as well k = = dan A = Z 
in the second term it is obtained, breaker length index of 


Hp _ 2tanh 6m 


F = ...(8) 
V. Equation of Breaker Depth Index “2 
It is defined as, 
A A 
tanhk (h +2) = pk (n+4) 
An equation is obtained as, 
a-p) =0 
By using (6) thus, 
_ tanh 67 
— Orn 
With this equation, the breaking equation becomes, 
Hp _ 4 tanh 67 9) 
hp  Ər-tanrhðr —s ( 
” hp 
VI. Equation for A 
hp 


The equation for F is formulated by using the 
b 


conservation equation of the wave number. 

In solving the Laplace equation with the variable 
separation method, the velocity potential is considered to 
be a multiplication of 3 functions (Dean (1991)), that is: 
pz t) = XTE) 

where X(x) is only the function of x, Z(z) is only the 
function of z and T(t) is only for the t time. In (1), which 
is Z(z): 

Z(z) = coshk(h + z) 

where h is for water depth. As the function is only of z 
then at sloping bottom where water depth h = h(x), also 
fork = k(x), 


dZ(z) ; dk(h + z) 
“ae T sinh k(h + z) ————— = 0 


dx 
In this equation the value is zero: 
dk(h +z) 
dx oe 
If this equation is worked on z = thus, 
ak (n +4) E 
dx 
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Thus, 
A 
k (n + 5) = constant 


By using (6), the equation for the conservation of the 
waave number is obtained, 


k(h+ 5) =9 
z) = OF 
The equation for the conservation of the wave number 


applies to all water depths, including the breaker depth, 
therefore it applies to the breaking point, 


Ap 
kp (hs + >) = On 


Substitute k, = 2” and Ap = 2b ' 
Lp FA 
L, 2 4L, 
Substitute = with (8), 
b 
Bee (10) 
Lp 2 21 
VII. | CONSISTENCY CHECKING ~, DAN we 
b b b 


The consistency test is proof that a breaker index is the 
product of the multiplication of the other two breaker 
indexes. The consistency test is done using the 
connectivity equation, 

2 = A uae (11) 
This equation states that the multiplication between the 


. H h ; H : 
equations = and - must create an equation of i which 
b b b 


is the same as (8). 


Substitute He with (9) and Pp with (10), 
hp Lp 


H, ( 4 tanh 0r )6 mn) 


Ly Om — tanh Or) \2 21 
It is obtained, 

H, 2tanh 8r 

Ip mtm 


The equation is the same as (8). Therefore, the equations 


Hp H h i : ; 
of ra and = meet the consistency requirements. This 
b Lp b 


consistency character also shows that the value of one 
breaker index is determined by another breaker index, or in 
other words there is interdependence between breaker 
indexes. Therefore, the formulation of the three breaker 
indexes should be done simultaneously so that consistency 
can be checked and there is connectivity between the 
breaker indexes. 


VIII. Breaker Height Index a. 
0 


The wave energy at one wavelength for a sinusoidal wave 
is 
E = cepgH*L (m) 
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: pa ea a 1 
Cg is a coefficient, which in linear wave theory cg = = 


(Dean (1991)). p is water mass density, g is gravitational 
force, H is wave height and L is wavelength. 


Based on the law of conservation of energy, the wave 
energy at the breaker point is the same as the wave energy 
in deep water, 

cepgHi by = crpgHĝ Lo 

The same elements cancel each other out, 

HÈL, = Bela 

Substituting the breaker length index (8), we get the 
breaker height equation associated with the breaker length 


index value 
_ 2tanh 67 


Hy = —— Hg Lo 


TT 


...(12) 
Hutahaean (2022) obtained the deep water wave number 
is, 


ko = Y3 


H 
(+B) 
y2 and y3 are coeficients where y, = 1.4 and y3 = 1.8. 
By this wave number, the deep water wavelength is: 


Lo = n(y2+¥2)Ho 
Y3 
Substituting Lo to (12), it is obtained 
1 
Hp 2 tanh(67)(y2 +13) /3 
==. << ..(13) 
Ho Y3 
IX. DETERMINATION OF VALUE OF 0 
The calculation results for Hb Hb he and Hp for some deep 
hp Lp Lp Ho 


water depth coefficient values 0 is presented in Table (2) 
below. To save space, only calculation results are 


; Hp . : 
presented where there is a match between a with previous 
b 


research. 


Table (2) the value 0 which results Ze that is in a line with 
b 


previous research 


0 H, Hp hp H, 


1.95 0.637 0.78 0.816 1.367 


2.1 0.637 0.715 0.891 1.367 


2.6 0.637 0.56 1.141 1.367 
There is a match between the values of ae from several 
b 


previous studies, where the match can be found in the 
value of @ that varies. 


-In @ = 1.95, it is obtained 2 = 0.78, which is in a line 
b 
with: 


a. Mc Cowan (1894) : z = 0.78 
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b. Weggel (1972) : 


21.56 
Hb _ IT /[1+exp(-19.5m)] 
hp 9 T2+hp43.75[1-exp(-19m)] 
i H 
m is bottom slope, for m = 0 Ta = 0.78 
b 


- In 8 = 2.10, it is obtained Ze = 0.715, which is in a line 
b 


with: 
a. Galvin (1969) : a ; for bottom slope 
hp 14-6.85m 
m < 0.07 
Inm=0, = + = 0.714 
hp 1.4 


b. Collins and Weir (1969) : =e = 0.72 + 5.6m 
b 


Inm = 0, 2 = 0.72 
hp 


c. Madsen (1976) : aad = 0.72(1 + 6.4m) 
b 


Inm = 0, Œ = 0.72 
hp 
- In 8 = 2.60, it is obtained a = 0.56, that matches with 
b 


Smith and Krauss (1989) 
H > — 1? 5(1 — exp (—43m)) = 
(0) 


hp = 1+exp(-—60m) 


H 
Inm = 0, = 0.56 
hp 
: hp . F 
To determine the @ the value of > is examined where 
b 


h : 
= <1 shows that the waterdepth is smaller than 
b 


wavelength. This condition causes water depth unable to 
support wave hydrodynamics resulting in breaking. In 0 = 


2.6, the value of ve > 1, it can be assumed that breaking 
b 


will not occur in that condition. 


Then, to select between @ = 2.1 and 0 = 1.95 where both 


. h : i 
result in = < 1, a review will be carried out based on deep 
b 


water depth criteria by using the wave number 
conservation equation in deep waters. Equation (6) is done 


in deep water, 
ko (ho + “2) = On 


2 


; 2 : 
Substitute k = and the equation can be expressed as 


Lo 2 4o 

ho in this case is the transitional depth between shallow 
water and deep water, where the water depth is smaller 
than hp is shallow water while the water depth is greater 
than họ belongs to deep water. At the depth of the 


ne h i 
transition, T should be greater than one or at least 1. This 
0 


suggests that 0 must be greater than 2. From this, it can be 
concluded that the breaker index that corresponds to Table 
(2) is the breaker index formulated by 0 = 2.1. 
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Breaker depth index in 0 = 2.1 produces a breaker depth 
index that is in accordance with the results of research 
from Galvin (1969), Collins and Weir (1969) and Madsen 
(1976). In addition, the values in Table (1) show that 
tanh 67 untuk 0 = 2.1 has reached 1 in the accuracy of 5 
decimal places while in 0 = 1.95 it has not reached 1. 
Both of these strengthen the conclusion that the 
appropriate value of the depth coefficient 8 is 2.1, with the 
breaker index on Table (2). 

From the above discussion it can be seen that the breaker 
index is determined by the deep water coefficient @ which 
means that the breaker index is set by determining the deep 
water depth. Thus there is a possibility that the diversity of 
breaker index obtained from previous research is also 
caused by differences in determining deep water depth. 

It has been shown that there is a match between the 
breaker depth index obtained with the breaker depth index 
from previous studies. In the following section, the value 
of the breaker length index will be studied = where the 


H . ; . H 
value of = obtained is a constant for all 0 that is r = 
b b 


0.637. This value is much larger than the critical wave 
steepness: 


a. Michell, J.H. (1894) : Z = 0.142 
b. Toffoli et al (2010) : z = 0.170 


Then it is studied the critical wave steepness from Michell, 
J.H. and Toffoli et al for the associated breaker depth 


index values. By specifying the value of = then @ can be 
calculate with (8), with the value of 0 obtained - using 
(9). The result for = = 0.142 — 0.170 is shown in Table 
(3). 


Table (3) the value 0 in some values of 


H — 0.142 — 0.170 
Lp 


Hp 0 Hp tanh 07 
Ly hy On 

0.142 0.072 | 233.951 | 0.98319 
0.146 0.074 | 220.914 | 0.98222 
0.150 0.076 | 208.905 | 0.98121 
0.154 0.079 | 197.819 | 0.98018 
0.158 0.081 | 187.565 | 0.97912 
0.162 0.083 | 178.060 | 0.97803 
0.166 0.085 | 169.234 | 0.97691 
0.170 0.087 | 161.023 | 0.97576 


In ab from 0.142 to 0.170, it is obtained the value of @ 
b 


which is smal, as much as 0.072-0.087, with =— ranges 
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from 0.98319 to 0.97576 closer to 1. The value of a 


obtained is very large not found in the results of previous 
studies. The large value is due to tanh 0m ~ Or. 


According to (6), in small amplitude it applies Om = kh 


tanhkh 


thus = 1, this is a shallow water long wave 


condition. From this condition it can be concluded that 
both Michell's (1894) and Toffoli et al's (2010) criteria are 
for long waves in shallow waters. 


a. Miche’s Equation Study. 
Miche (1944), create a breaker index equation connecting 


the two breaker indexes, that is 
H 


b_ 2mhp 
T = 0.142 tanh (= ee (14) 
The equation connects two breaker indexes of = and a 

b b 


ta) H h 

By determining the value of oa thus the value of = can be 
b b 

calculated. After that, by using the connectivity equation, 


H . : : 
= can be calculated. The calculation results is shown in 
b 


Table (4). 
Table (4) Calculation of the breaker index using the Miche 


equation. 
Hp hp Hp 
Ly Ly hy 
0.082 | 0.10483 0.78224 
0.099 | 0.13716 0.72179 


0.126 | 0.22428 0.5618 
0,142 | 1.15351 0.12397 


; . H 
In Table (4) it can be seen that there is a match between 4 
b 


; ; . : Ses 
and previous studies and in = maximum which is 0.142. 
b 


; n H h 
It was obtained that the match is in the values of = and = 

b b 
is very small, in which this is the condition of long wave in 
shallow water. 


Then, the Miche’s equation is modified to, 
= Z 0.637 tanh (=) 
L 


b Lp 
The calculation results where there are matches between 


(15) 


Z and previous studies are presented in Table (5). 
b 


Table (5) Calculation of breaker index using (15). 


Hp hp Hp 

Ly by hy 
0.636955 | 0.815732 0.780839 
0.636980 | 0.880247 0.723638 
0.636995 | 0.990467 0.643126 
0.637000 | 1.153513 0.552226 
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Table (5) indicates that there is conformity with Table (2) 


; Hp h H ; i 
in the value of a - and also = . Then Miche's equation 
b Lp b 


can be modified into an equation for short waves by using 

H : H : ; 

= 0.637 or in other words, aa = 0.637 is according to 
b b 


Miche's equation. 


X. CALCULATION RESULTS OF BREAKING 
PARAMETERS USING THE BREAKER INDEX 


H H 


In section 9, it was obtained = 0.637, = 0.715, 
b b 

æ = 0.891 and = = 1.367 in 0 = 2.1. An example of 

b 0 


calculating the breaking parameter with the breaker 
indexes for several deep water wave heights Hp) is 
presented in Table (6). 

Table (6) The calculation of parameters breaking 


Ho T H, Lp hy 
(m) (sec) (m) (m) (m) 
1 4.61 1.37 | 2.15 1.91 


1.5 5.65 2.05 3.22 2.87 
2 6.53 2.73 4.3 3.83 
2.3 7.3 3.42 5.37 4.78 
3 7.99 4.1 6.44 5.74 


Wave period in Table (6) is obtained from the relation of 
wave period with wave height (Hutahaean (2022)), 


872 (v2 +> 


2a (16) 


v3)" Ho 


T= 


To get an overview of the conditions resulting from the 
calculation of the breaker parameter with the breaker index 
equation obtained, we will review the breaker height from 
the previous breaker height equations. 


The empirical breaker height index equation F is quite a 
0 


lot and can be divided into two groups, namely those that 
use the bottom slope as a parameter and those that do not 
use. In this study, the breaker depth index is used without a 
bottom slope as a parameter and or the bottom slope is 
given a value of zero. Researchers in this group include 
Komar and Gaughan (1972), Larson and Krauss (1989), 
Smith and Krauss (1990) and Gourlay (1992). 

Komar and Gaughan (1972) : 


Hp _ Ho\ 722 
ab = 0.56 (=) ae (17) 
Larson and Kraus (1989 : 
Hp _ Ho) 70-24 
ab = 0.53 (=) ee (18) 
Smith and Kraus (1990) : 
Hp _ Ho\70:3040.88 m 
ab = (0.34 + 2.47 m) (=) (19) 
m is the bottom slope, this research uses m = 0 
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Hy Ho\~928 
Gourlay (1992) : $Ë = 0.478 (=) (20) 
From these equations, the deep water wave length is 


calculated using the dispersion equations of linear wave 
_ 9 

theory, Ly = = 

The breaker height calculation results using the empirical 

breaker height index equations are presented in Table (6). 


Table (7) Breaker height from empirical equation. 
Hy(m) 


Ho T 
(m) | (sec) | (17) (18) (19) (20) 
1 | 461 | 113 | 1.23 | 0.97 | 1.27 
15 | 5.65 | 169 | 184 | 146 | 1.91 
2 | 653 | 2.26 | 246 | 1.95 | 2.55 
2.5 | 73 | 282 | 3.07 | 243 | 319 


3 7.99 3.39 3.69 2.92 3.82 
Note : (17), (18), (19), (20) are the equation code number. 


Of the four breaker height index equations, (20) gives the 
largest and closest result to the breaker height in this study. 
Furthermore, the difference between the breaker height 
from this study and the breaker height from (20) was 
studied. The comparison results are presented in Table (8). 
Table (8) Comparison between breaker height and 
Gourlay’s breaker height 
Ho T Hp (m) 5 
(m) | (sec) | (3) | 20 | (%) 
1 4.61 1.37 | 1.27 7.24 
1.5 5.65 2.05 | 1.91 7.24 
2 6.53 2.73 | 2.55 7.24 
2.5 73 3.42 | 3.19 7.24 
3 7.99 4.1 3.82 7.24 


Note : 5 = | x100 % 


In Table (8), the breaker height of the results of this study 
is greater than the breaker height of (20) with a constant 
difference of 7.24%. To bring the results (13) closer to 
(20), it can be done by reducing the value of 8, when 0 = 


0.35, it is obtained ô = 0.44 %, but if ae = 10.7, it shows 
b 


that the (20) is for long wave. 


XI. CONCLUSION 
The first conclusion from this study is that the Kinematic 
Free Surface Boundary Condition has breaking 
characteristics and the breaker index equation can be 
formulated. Meanwhile, in the wave number conservation 
equation there is a breaker index in the form of a ratio 
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between breaker depth and breaker length. The breaker 
height equation can be obtained by substituting the breaker 
length index in the energy conservation equation. With this 
procedure the breaker index equations for short waves are 
obtained. 

From the results of the connectivity test on the resulting 
breaker index equations, it was found that there is a link 
between the breaker indexes. This also gives the 
conclusion that there is interaction between breaking 
parameters that the value of one breaking parameter 
depends on the value of another breaking parameter. 

Thus it can also be concluded that breaking wave research, 
both laboratory and analytical research cannot be done 
separately, it should be done simultaneously for all breaker 
indexes or all breaking parameters. 

Determining the deep water depth affects the breaking 
parameters, especially the breaker depth. Therefore it is 
estimated that the difference between the breaker depth 
index results of previous studies is due to differences in the 
determination of deep water depth. 
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